Role of vibrations in the jamming and unjamming of grains discharging from a silo 
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We present experimental results of the jamming of non-cohesive particles discharged from a flat 
bottomed silo subjected to vertical vibration. When the exit orifice is only a few grain diameter 
wide, the flow can be arrested due to the formation of blocking arches. Hence, an external excitation 
is needed to resume the flow. The use of a continuous gentle vibration is a usual technique to ease 
the flow in such situations. Even though jamming is less frequent, it is still an issue in vibrated silos. 
There are, in principle, two possible mechanisms through which vibrations may facilitate the flow: 
(i) a decrease in the probability of the formation of blocking arches, and (ii) the breakage of blocking 
arches once they have been formed. By measuring the time intervals inside an avalanche during 
which no particles flow through the outlet, we are able to estimate the probability of breaking a 
blocking arch by vibrations. The result agrees with the prediction of a bivariate probabilistic model 
in which the formation of blocking arches is equally probable in vibrated and non-vibrated silos. 
This indicates that the second aforementioned mechanism is the main responsible for improving the 
flowability in gently vibrated silos. 

PACS numbers: 45.70.-n 



I. INTRODUCTION 



Granular materials present interesting and unusual 
physical properties [3, [2| ■ As the collisions among grains 
are highly dissipative, dynamical states (such as a gran- 
ular flow) need a continuous energy supply. Otherwise, 
dissipation quickly stops the motion. Even with a sus- 
tained energy input, granular flows can -in some cases- 
get jammed and halt [H, H|. 

A classical example of a system prone to jamming is a 
silo discharged by gravity. The flow of particles may be 
suddenly arrested by the formation of a blocking arch (or 
vault) if the size of the outlet does not exceed a few times 
the size of the particles 0, S 0, S S M H3 • 

Controversy 

still goes on about whether there exists or not a certain 
size of the orifice above which jamming is not possible. 
It seems that such a critical size can be found for a 3D 
silo but not in a quasi-2D setup [l(| EH ■ Whatever 
it may be, in many cases it is interesting or unavoid- 
able to reduce the size of the orifice in order to limit 
the flow rate; and the price to pay is an increased jam- 
ming probability. Hence the significance of the methods 
used to improve the flowing of grains, among which vibra- 
tions (either local or global) applied to the silo are quite 
common [H, [H [3 M, M, Ed, E 03 ■ Unfortunately, 
clogging still occurs. During the discharge of a vibrated 
silo through a small orifice, blocking arches may form. 
These structures can stop the flow even when the silo is 
being vibrated. However, the flow may be spontaneously 
restarted thanks to the destabilization of a blocking arch 
induced by the continuous excitation. Eventually, some 
robust jam will develop, in the sense that the continuous 
vibration is unable to restart the flow. A stronger exter- 



nal perturbation is then needed to unjam the system. 

Previous works [1, in non vibrated silos character- 
ized the jamming by measuring the amount of grains 
delivered between two consecutive jams. For lack of a 
better term to designate this quantity, it was called an 
avalanche. We have borrowed this word because it con- 
veys the idea of a granular flow that starts and stops 
abruptly, but should not be confused with the rapid land- 
slide or flood of material along a slope, and in particular, 
on the surface of granular piles [2CJ, |2l| . In the discharge 
of a silo, the avalanche size is the number of grains that 
flows through the orifice from the moment when the out- 
pouring starts until the formation of a blockage that ar- 
rests the flow. In the absence of vibration, the avalanche 
size is a well defined quantity, as the arches or vaults that 
block the orifice are robust: once the flow is stopped, it 
does not resume by itself even after a long waiting time. 
It has been shown that in silos of two and three dimen- 
sions the avalanche size distribution displays an expo- 
nential tail. This can be understood if each grain has a 
probability p of passing through the outlet without form- 
ing a blocking arch that is constant during the avalanche 
@, H|. In a recent paper Janda et al. [ll[ related the 
value of this probability in a quasi-2D silo to the arch 
size distribution within a 2D static granular layer. 

When vibrations are applied, the notion of avalanche is 
further complicated due to the fact that arches blocking 
the orifice can break down after a certain time. This re- 
sults in an intermittent flow; in this case, in order to 
conclude that the flow has stopped permanently, one 
has to wait for a long time after grains cease to come 
out from the silo. The avalanches can contain them- 
selves short lapses during which the flow has temporarily 
stopped, due maybe to arches that blocked the exit and 
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FIG. 1: Experimental device. S: silo; V: pneumatic vibrator; 
N: damping and isolation system (air cushions and valves); P: 
photosensor; A: compressed air duct; B: box; S: scales. The 
inset shows a photograph (from the bottom) and a sketch of 
the cross section of the changeable part in which the hole is 
bored. 



were broken apart by the vibrations shortly afterwards. 
In Sect. IHII we will explain how we choose a meaningful 
waiting time for the jams. 

In this work, we present experimental results on the 
avalanche size distribution that allows for a comparison 
between vibrated and nonvibrated silos. We are able to 
provide the jamming probability in vibrated silos, and 
formulate a simple model that extends previous results 
and allows us to quantify the effect of vibrations. It can 
be conjectured that the application of vibrations to a 
silo may lead to two mechanisms that could cause a re- 
duction of robust jamming events: (i) a decrease in the 
probability of the arch formation, which are the cause 
of blockages, and (ii) the loss of stability of these arches 
after they are generated. We will show that the second 
mechanism is the responsible for the appearance of larger 
avalanches when a gentle vibration is applied. 



II. EXPERIMENTAL SETUP 

The experimental setup (Fig. [I} comprises a flat bot- 
tomed cylindrical silo with a circular hole of diameter D 
in its base. The part making the outlet orifice can be 
replaced, and in this way D can be changed. The silo is 
filled with Delrin spheres (diameter d = 3.00 ± 0.01 mm 
and weight w = 18.9±0.03 mg). The dimensionless size 
R of the orifice is defined as R = D/d. The diameter of 
the silo is more than 30 times the particle diameter in 
order to prevent that the lateral walls influence the jam- 
ming probability [22j ■ The silo is routinely refilled so that 
the height of the granular material is always more than 
twice the diameter of the silo. A pass-through photosen- 
sor was placed just below the outlet in order to detect 
whether beads are flowing or not at every time; the res- 
olution is such that one single bead cannot go through 
it undetected. At the bottom of the silo, a box collects 
the grains that fall from the exit. This box is placed on 
an electronic scales, with an accuracy such that a sin- 
gle bead can be resolved. When the flow is arrested (see 
below the criterion we use to define that the flow has 
stopped permanently), the weight of the avalanche is ob- 
tained from the scales, and its size in number of beads is 
calculated. Then, a new avalanche is triggered with a jet 
of compressed air from beneath the orifice, a technique 
used in previous works [1, Q . 

The whole silo is continuously shaken by a pneumatic 
vibrator actuating in the vertical direction. In order 
to facilitate the vibrating movement and to isolate it 
from the measuring devices, the structure is supported by 
three air cushions, controlled through mechanical valves. 
We have carried out the measurements at a fixed effec- 
tive dimensionless acceleration T = a e ff/g — 0.22 and at 
a frequency of 110 Hz. Here a e // is the rms value of the 
acceleration as measured by an accelerometer attached 
to the silo base, averaged over 200 cycles. In fact, this 
experimental set-up is quite similar to another used in 
previous works (see [9]) with the addition of the pneu- 
matic vibrator and the air cushions on which the silo 
rests. 



III. AVALANCHE SIZE DISTRIBUTION AND 
JAMMING PROBABILITY 

It is interesting to begin by inspecting a typical signal 
from the photosensor, which is binarized to one or zero 
depending on whether or not a particle is blocking the 
light beam. Two signals (each one registered over the du- 
ration of a single avalanche) are displayed in Fig. [2] one 
of them was obtained when the silo was being vibrated 
(Fig. Ob), and the other was taken with the vibration off 
(Fig.[5]a). In the latter case, there are short interruptions 
of the flow; but once the orifice gets jammed, the flow is 
arrested permanently. On the other hand, when the silo 
is being vibrated, there are stretches during which grains 
are flowing separated by time intervals when the flow is 
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FIG. 2: Signal from the photosensor at the silo exit: a value 
of one indicates that a particle is blocking the beam, zero 
means that the beam is unobstructed, (a) Nonvibrated silo, 
(b) Vibrated silo, (c) A zoom of the signal shown in (b) 
during the first three seconds, the same time stretch than in 
(a). All the data shown have been obtained for an orifice size 
i? = 3.05. 



interrupted. These can be significantly longer than in 
the nonvibrated case. This means that a jam has devel- 
oped, but the vibrations are able to restart the outpour- 
ing. Eventually a blockage is formed that stops the flow 
permanently because the external vibration cannot break 
down the blocking arch; at least within our experimental 
time scales. In this sense, Fig.[5]b could be described as 
if there were avalanches inside the avalanche. Fig. [2] c is 
a zoom of Fig.[2]b that reveals that the form of the signal 
of the vibrated silo, in a short interval during which the 
material is flowing, is similar to the one displayed by the 
nonvibrated silo. 

As we mentioned in Sect. HI the avalanche size s is the 
number of particles fallen between two consecutive jams. 
For the nonvibrated silo, the end of the avalanche is eas- 
ily detected, as once the silo is jammed the flow does 
not restart by itself. In the vibrated case, we consider 
that the silo is jammed whenever the outflow stops for a 
time longer than 100 s. We have seen that is extremely 
rare, for the vibrations used in our experiments, that the 
silo resumes its discharge after being jammed for more 
than this amount of time. Up to a certain point, this 
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FIG. 3: Probability n s (R, F) of finding an avalanche of size 
s for an orifice size R = 3.02, without vibration (r = 0, ■) 
and with vibration (F = 0.22, o). The solid lines correspond 
to fits using Eq. ([4]). For the nonvibrated silo, the fitting 
parameters are p = 0.981 and q — 0. For the vibrated silo 
p = 0.981 and q = 0.836. 



value is arbitrary (in the sense that one could have de- 
cided to wait for a longer time). However, the choice was 
made because the results presented here do not change 
significantly if a sufficiently long time cutoff is taken. Of 
course, it must be understood that this value can change 
a little with experimental conditions, such as T. 

For each value of R, the size of about 3000 avalanches 
has been measured for the vibrated and the nonvibrated 
silo. We define n s (R,T) as the probability of finding an 
avalanche of s grains in a vibrated silo with an opening of 
size R and acceleration T. As can be seen in Fig. [3J the 
PDF of the avalanche size displays an exponential tail for 
large avalanches. The comparison between the PDFs for 
the vibrated and nonvibrated cases (shown in the same 
figure for a particular value of R) reveals that the use 
of vibrations increases the size of the avalanches. How- 
ever, the shape of the distribution remains exponential 
for large avalanches. We have found this behavior for all 
the outlet sizes explored (1.7 < R < 4.0). We emphasize 
here that jamming has not been observed for R > 5 in 
the nonvibrated silo @. 

The jamming probability J/v(i?,r) is defined as the 
probability that the silo jams before N particles fall 
through the orifice for given R and T. In other words, 
Jn(R,T) is the probability of finding an avalanche 

smaller than N, i.e. J N (R,T) = J2s=o n ^ T ) 0- In 
Fig. 2] we show the instance Jioo(R, T), obtained from 
the experiments with both the vibrated and the nonvi- 
brated silo. Jioo(-R, r = 0) is close to one for small val- 
ues of R and falls to zero rather sharply in the interval 
2.8 < R < 4.0. In 3D silos it has been shown that for 
N — > oo the fall tends to a Heaviside function centered 
in i?« 5 [1,0. 
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FIG. 4: Jamming probability Jioo(R, F), i.e. the probability 
that the silo gets jammed before 100 particles flow through 
the outlet, as a function of R, for the vibrated (o) and the non- 
vibrated (■) silo. The lines correspond to a fit using Eq. J2). 
The "weak jamming" probability Jxoo(R) f° r the avalanches 
within an avalanche is also shown (A). 




R 

FIG. 5: Mean avalanche (s) versus the dimensionless radius R 
for the vibrated (o) and nonvibrated (■) silo. The solid lines 
correspond to fits using Eq. (TJJ with the parameters given in 
the text. The dotted line and the dashed line correspond to 
the values of R c used in the fits. 



IV. PROBABILISTIC MODEL FOR ARCH 
FORMATION AND BREAKAGE 



The vibrated case displays the same trend; however, 
the sharp decrease in J 100 (R,T = 0.22) occurs at lower 
values of R (2.5 < R < 3.5). It is clear that the vibration 
particularly reduces the jamming probability in the range 
of opening sizes where the transition from high to low 
jamming probability occurs. 

For a stationary silo, the functional dependence of the 
mean avalanche size (s) with R was shown to be compat- 
ible with a power law divergence Q : 



C 



(Rc - RY 



(1) 



Let us assume this dependence also for vibrated silos. 
The experimental data can be fitted with this equation, 
as shown in Fig. [5] The values obtained for the fitting 
parameters are C = 2.3 ± 0.7 x 10 5 , R c = 5.52 ± 0.03 
and 7 = 8.62 ± 0.19 for the nonvibrated case, whereas 
C = 1.21 ± 0.3 x 10 4 , R c = 4.6 ± 0.1 and 7 = 7.56 ± 0.52 
for the vibrated case. 

The functional form of the jamming probability 
Jpf(R,T) can be related to (s). This was done for a non 
vibrated silo, as explained in [llj. Hence, using Eq. |T|), 
an expression of J/v(-R,r) can be obtained: 



Jn(R,T) 



-N(R C -R)~< /C 



(2) 



which reveals an excellent agreement with the experimen- 
tal data (see Fig. [3]) if the values of the fitting parameters 
obtained from Eq. (|TJ) are used. 



In order to describe the intermittent flow in a vibrated 
silo, one can assume that arches may form at the exit just 
as if the vibrations were absent. If the vibration breaks 
down the arch, the discharge resumes. Eventually, an 
arch strong enough can appear such that the vibration is 
not able to break it, at least for a considerable amount of 
time. As explained above, if an arch lasts for more than 
100 seconds, it is assumed that it is robust enough and 
the flow is then said to have been stopped. Alternatively, 
it could also be conjectured that somehow the vibration 
inhibits the formation of arches. A reduced probability of 
arch formation would also result in longer avalanches. We 
now focus on answering the question of whether a gentle 
vibration diminishes the arch formation probability or, 
instead, operates by destabilizing blocking arches formed 
at the exit. Of course, a mix of the two mechanisms 
may be in place; we can still investigate if one of them 
dominates. 

Let us extend a previous model Q to the case of a 
vibrated silo. We introduce a bivariate probability dis- 
tribution to represent the behavior of a vibrated silo. We 
define p(R, T) as the probability that a single grain passes 
through the outlet without forming a blocking arch along 
with its neighbors during continuous flow. (In the follow- 
ing we will drop the functional dependence of p in the no- 
tation). Therefore l—p accounts for the probability that 
a particle does get involved in the formation of an arch 
as it reaches the outlet. This probability p is assumed to 
be the same for all grains and independent of discharge 
history. 

Now, we introduce a new variable q(R, T), which repre- 
sents the probability that, once a blocking arch is formed, 
a particle flows through the exit as a consequence of the 



5 



destabilization of the arch due to the continuous vibra- 
tion. The probability q is a measure of the "arch break- 
age probability" , and does not take into account the time 
that destabilization may require for a given arch. There- 
fore, the probability that a blocking arch remains stable 
— i.e., no more particles will ever flow out of the silo — 
will be 1 — q. We assume that all arches blocking a given 
orifice have the same breakage probability q for a given 



J 



vibration intensity. This is a fair approximation since, 
for a particular R, the blocking arches are expected to 
have a typical size. Notice that this might not be true in 
hoppers, where arches form at different positions along 
the hopper walls, and hence a wide distribution of arch 
sizes may be present for a given exit size. 

The probability n s (R, T) that an avalanche consists of 
s grains can be obtained as 



n s (R,T)=p s (l-p)(l-q) 



(s-0)!0! 



+ p s - 1 (l-p)q(l-p)(l-q) 



(s-l)!l! 



+ p s - 2 (l-p)q(l-p)q(l-p)(l~q) 



+ ... 



+ p s - k [(l-p)q] k (l-p)(l-q) 



(s-2)!2! 



(s - k)\k\ 



+ {(l-p)qY(l-p)(l-q) 



0!s! 



< — a single blocking arch forms and remains stable. 

< — an arch forms and breaks, then a new arch forms 
and remains stable. 

< — two arches form and break, an arch forms and 
remains stable. 

< — k arches form and break, an arch forms and 
remains stable. 

< — s arches form and break, an arch forms and 
remains stable. 

(3) 



In Eq. |3]) the factor sl/((s — k)\k\ in each term ac- 
counts for the number of permutations with repetitions 
of the s — k events of probability p (in which a grain 
flows through the orifice when there is no blocking arch 
present) and the k events of probability (1— p)q (in which 
a grain forms a blocking arch with its neighbors followed 
by the flow of a grain due to arch destabilization) . These 
events can happen in any order but always add up to 
the desired number s of grains flown. Each of these per- 
mutations will correspond to a set of k jams that are 
unstable under the applied vibration, separated by the 
continuous flow of grains. For example, the last term in 
Eq. ([3]) corresponds to the extreme case where the flow 
of every single particle is followed by the formation of an 
arch that then breaks due to the vibration. As mentioned 
above, the avalanche finishes when the last clogging hap- 
pens and remains stable [probability (1— p){l — q)]. From 
Eq. ©, n s (R, T) can be written as: 



n s (R, r) is a decreasing exponential function of s, in ex- 
cellent agreement with the experimental results shown in 
Fig. [31 We can see that n s is properly normalized, i.e., 
^^ n s (i?) = 1, since this sum is a geometrical series 
with ratio r < 1 and it converges to one. 

Thus, we can calculate the mean avalanche size, de- 
fined as 



< s > (R, r) = J2 sn s(R)- 



(5) 



s=0 



By inserting Eq. (01 in Eq. ([5]) we obtain a Gabriel's 
staircase [23. Therefore 



< s > (r, r) = 



p+ (i -p)q 
(i-p)(l-g) 



(6) 



,(R,r) = (i-p)(i-q)J2p s - k [(i-p)q] k ] 



k=0 



(s - k)\k\ 

= (l-p)(l-q)[p+(l-p)q} s . (4) 
It is simple to show that p + (1 — p)q < 1. Therefore, 



V. ESTIMATE OF THE ARCH BREAKAGE 
PROBABILITY 

Equations (0| and ^ reduce to the expression for the 
nonvibrated case [ll| setting q = 0. Equation Q yields 



0.0 -I 1 ' 1 1 1 ■ , ■ 1 ' r- 

1.5 2.0 2.5 3.0 3.5 4.0 

R 



0.01 



Ol ^~ 1 
At (S) 



10 



FIG. 6: Probability q that a blocking arch is destabilized 
due to the vertical vibration as a function of the opening size 
7?. Results obtained with the probabilistic model, in which p 
does not depend on F (solid circles). Results obtained from 
the analysis of the weak jams (open squares). The data cor- 
respond to several sets of experiments; for some of them both 
methods were implemented, while for others only one of the 
measurements was performed. 



FIG. 7: Normalized histograms of the time intervals (At) 
within an avalanche during which there are no particles flow- 
ing through the orifice. Closed (open) symbols represent the 
data obtained with the nonvibrated (vibrated) silo for the val- 
ues of R indicated in the legend. For each R, measurements 
have been carried out within, at least, 3000 avalanches. 

open squares in Fig. [5] correspond to an independent 
estimate of q measured as explained below. 



an exponential form of n s for both vibrated (q ^ 0) and 
nonvibrated (q = 0) silos. In Fig. [3] we show the fit of 
the results for the nonvibrated case using Eq. g]) with 
q = and p = 0.981. For the vibrated case the data 
are fitted using Eq. ^ with the same value of p but 
q = 0.836. This suggests that p may well be independent 
of r, at least for gentle vibration, which in turn implies 
that the probability 1 — p of arch formation is not affected 
by vibration. Hence, the enlarged sizes of the avalanches 
obtained for the vibrated silo would be due solely to the 
appearance of a probability q associated to arch destabi- 
lization. 

In what follows we estimate q assuming that p does 
not depend on the external vibration. According to Eq. 
©, we can obtain p(R, T = 0) from the mean avalanche 
size of the nonvibrated silo as p(R) = p(R, T = 0) = 
( s )_R.r=o/[l + ( s )i?,r=o]- Thus, from Eq. ([6]), q can be 
calculated for T > with the mean avalanche size for 
both cases as: 



q(R,T^Q) = -f — (7) 

In Fig. [6] we represent with solid circles the values of 
q obtained through Eq. ([7]) by using the measured val- 
ues of (s), for a vibrated and a nonvibrated silo. It is 
not surprising that q increases with i?, since big arches 
-needed to block big orifices- are expected to be less ro- 
bust against small perturbations than small arches. The 



VI. THE INTERNAL STRUCTURE OF 
AVALANCHES 

It is possible to obtain an additional, independent esti- 
mate of the arch breaking probability, motivated by the 
observation of a behavior that has already been pointed 
out (Fig. [5]): during the flow of an avalanche, short clogs 
are registered. We have measured the time intervals At 
during which no beads cross the exit orifice. We repre- 
sent in Fig. [7] the normalized distribution of At for two 
different values of R, and these for a vibrated and a nonvi- 
brated experiment. Let us remark that it is not judicious 
to consider intervals much shorter than the time it takes 
for a bead to traverse the orifice due to gravity (which is 
t = \J Wl '<?? about 0.02 s for the grains considered here). 

The first noticeable feature in Fig. [7] is that for a non- 
vibrated silo, the probability of finding a time interval 
during which the flow is temporarily stopped decreases 
abruptly with At (note the logarithmic scale in the plot). 
With vibrations, things are different. The histograms of 
Fig. [7] display a finite number of events for large At, in- 
dicating that within the avalanches longer intervals are 
present during which no grain comes out from the exit 
orifice. These events were not observed in the nonvi- 
brated silo (nor in recent results obtained in a 2D silo 
[2~il]). These long events suggest that there are flow in- 
terruptions that might have resulted in the end of an 
avalanche should the silo had not been vibrated. The 
distribution of At smaller than 0.1s shows similar behav- 
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ior for vibrated and nonvibrated silos. Clearly, this part 
of the distribution corresponds to the temporal scales 
associated with the "interparticle" time intervals in the 
continuous flow [24[ • Figure suggests that a temporal 
cutoff can be established to separate the vibrated and 
the nonvibrated dynamics of a few tenths of a second 
(where the distributions for the vibrated and nonvibrated 
cases begin to differ). We decided to take At = 0.5s as 
this cutoff. Although this choice is somehow arbitrary, 
we checked that the results presented below do not de- 
pend strongly on this value (for example, At = 0.25s and 
At = Is yield similar outcomes). 

Following this line of reasoning, let us call a weak jam 
all the events corresponding to At > 0.5s. Then we 
can readily count the number of weak jams (J w ) and 
take that value as an indication of the number of times 
an arch has been formed and broken due to vibration. 
Since the number of robust jams (j r ) corresponds to the 
total number of avalanches studied for given R and T, 
q{R, r) can be obtained by dividing the number of weak 
jams j w by the number of total jams (weak plus robust): 
q(R,T) =j w /(j w + j T ). 

In Fig. E] the results for q(R,T = 0.22) obtained with 
this direct measure are reported (open squares). The 
agreement with the values of q obtained from the proba- 
bilistic model under the assumption of the independence 
of p from r is reasonable. We then claim that the prob- 
ability 1 — p that a particle forms a blocking arch along 
its neighbors at the time of passing through the exit dur- 
ing the continuous flow does not depend on the external 
vibrations applied to the system, at least for the gen- 
tle vibration intensity we used. Therefore, the main role 
of vibration is to break blocking arches once they are 
formed. 

It is now possible to measure the size of the avalanches 
inside the avalanches, with which we refer to the amount 
of material unloaded in the vibrated silo between two 
weak jams, as explained above. From these results 
we obtain the avalanche size distribution inside an 
avalanche, and the corresponding weak jamming proba- 
bility J^(R, r). Since we assumed that the cutoff chosen 
to identify weak jams is such that these truly correspond 
to the jams that would be stable if vibrations were absent, 
we expect the probability distribution of those avalanches 
in a vibrated silo to display the same properties as nor- 
mal avalanches in a nonvibrated silo. In Fig. Q] we can 
see this is indeed the case, since Ji 00 (R, T = 0.22) (A) 
coincides with J W0 (R,T = 0) (■). 

VII. CONCLUSIONS 

We have reported experimental results of the avalanche 
size distribution and jamming probability in the dis- 



charge of a vertically vibrated silo. The first result that 
becomes apparent is that the presence of vibrations sig- 
nificantly increases the size of the avalanches. Addition- 
ally, we have shown that in a vibrated silo the discharge 
of grains is intermittent, due to the appearance of tem- 
porary blockages of the exit orifice. Taking as a starting 
point a previous model for a nonvibrated silo, we have 
extended it and explained the new results by introducing 
a bivariate probability distribution taking into account 
(a) the probability p that a particle pass through the 
outlet without forming an arch with its neighbors during 
continuous flow, and (b) the probability q that a particle 
flows through the exit -once and arch has been formed- 
thanks to the destabilization promoted by vibration. As- 
suming that p does not depend on the presence of vibra- 
tions, which is a fair assumption for gentle vibrations, 
the values of the mean avalanche size for vibrated and 
nonvibrated silos provide a straightforward estimation of 
q for different values of R. Besides, a direct measure 
of q has been carried out by recording the time interval 
within each avalanche during which grains do not flow 
through the exit. The two measurements agree remark- 
ably well, indicating that the main effect of vibration is 
the breakage of arches once they are formed without af- 
fecting significantly the arch formation probability. 

Additionally, the result of the monotonic increase of q 
with respect to the outlet size displays an expected be- 
havior: the greater the arch, the more unstable against 
small perturbations. Further measurements of the sta- 
bility of arches must be carried out for different vibra- 
tion frequencies and strong amplitudes to asses the full 
range of vibration parameters for which this explanation 
is valid, i.e., a rigorous definition of what a gentle vibra- 
tion means must be provided. 
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